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Abstract 

We describe the statistics of chaotic wavefunctions near periodic orbits using a basis of 
states which optimise the effect of scarring. These states reflect the underlying structure 
of stable and unstable manifolds in phase space and provide a natural means of charac- 
terising scarring effects in individual wavefunctions as well as their collective statistical 
properties. In particular, these states may be used to find scarring in regions of the spec- 
trum normally associated with antiscarring and suggest a characterisation of templates for 
scarred wavefunctions which vary over the spectrum. The results are applied to quantum 
maps and billiard systems. 

1 Introduction 

The observation that random matrix theory (RMT) describes the statistical properties of the 
eigensolutions of a classically chaotic quantum system pQ has formed a pivotal role in our 
understanding of such systems. The most basic prediction for eigenfunctions is that they form 
a Gaussian distribution when projected on a generic state. More explicitly, the overlap between 
a generic probe state \ip) and the eigenstates \ip n ) of the chaotic system has the probability 
distribution function 

P(x) = -±= ^ 2/2 , (1) 
y An 

where we normalise the states \if> n ) so that the overlap x = (f\ip n ) has unit variance and assume 
time-reversal invariance for ease of notation. This distribution has been shown to work very 
accurately for billiard systems, for example [2|. 

In any such comparison with RMT statistics, deviations due to short time dynamics and 
other system-specific effects are particularly interesting. In systems such as quantum dots (HUE], 
for example, such deviations may help us distinguish between RMT statistics arising from a 
simple underlying chaotic dynamics and RMT statistics which arise simply because the system 
is complex, perhaps dominated by interactions etc. In the context of wavefunction statistics, 
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a remarkably strong example of such deviation has been demonstrated by Kaplan, Heller and 
coworkers El Ej- They find that wavefunction statistics around a periodic orbit of the 
classical limit can deviate strongly from (0) , an effect attributed to the phenomenon of scarring 
[El 13 [Jill EJ] (whereby wavefunctions of classically chaotic systems are seen frequently to behave 
atypically near periodic orbits). These deviations can be robust and survive the semiclassical 
limit if the probe state \(p) is sufficiently localised in phase space around a periodic orbit. In this 
paper we describe the collective statistics of many components of the wavefunction using a basis 
of probe states which are locally complete near a periodic orbit. We find that a particular set 
of probe states can be defined in terms of which these deviations may be described particularly 
simply and these are related to states which have previously been used to characterise scarring 
under different guises (0 El El and [TJ1 ECU EH EH EB1), known variously under the names of 
universal test states [Hj, scarfunctions [Hj or simply quasimodes [TB] . 

The first theoretical explanation for scarring [0] was based on the evolution of a Gaussian 
wavepacket launched on a classical periodic orbit. Fourier analysis of the evolving wavepacket 
shows that typical overlaps of the wavepacket with chaotic eigenstates ip n fluctuate quasiperi- 
odically as a function of energy. Peaks occur in scarred parts of the spectrum associated with 
enhanced wavefunction probabilities near the periodic orbit and troughs are found in the anti- 
scarred parts where wavefunctions are expected to be suppressed near the periodic orbit. More 
recently this basic idea has been developed further by Kaplan et al El Ej to produce a 
quantitative theory of wavefunction statistics near a periodic orbit. 

Kaplan et al consider the statistics of overlap with a single probe state \ip), which might 
typically be Gaussian. By considering quantum recurrence for such a state launched on a 
periodic orbit, one can show that the variance of the distribution p(x) must deviate from the 
constant value predicted by RMT and is a function of the spectral parameter — energy E in 
the case of time-independent systems and eigenangle < 6 < 2n in the case of quantum maps. 
The overlap distribution in (Q) is then naturally replaced by 



where we adopt the notation of eigenangle for the spectral parameter and where the variance 
C(8), which may be calculated from the linearised autocorrelation function of the probe state, 
is independent of Ti. 

This was generalised in to provide a joint-probability distribution describing the statis- 
tics of the wavefunction in a basis of probe states which are complete in a neighbourhood 
in phase space of the periodic orbit. A minimum- uncertainty Gaussian wavepacket can be 
thought of as the ground state of a harmonic-oscillator Hamiltonian. We can also use as probe 
states the excited states \<pk) of the same Hamiltonian and thereby obtain a series of over- 
laps x = (x ,Xi, ■ ■ ■ ,Xk, ■ ■ ■) representing the components of a chaotic eigenfunction in the 
harmonic-oscillator basis. (We adopt notation here appropriate to a single degree of freedom 
transverse to the periodic orbit, but it should be clear how to generalise these statements to 
higher dimensions.) One can then derive a matrix of correlations C(9) with components 





C k i(0) = (x k Xi) 
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which depends as before on the spectral parameter 9 and which is independent of K (and there- 
fore survives the semiclassical limit). The corresponding normal joint-probability distribution 
for the overlap amplitudes 

P(x;0) = m e-^"'^/ 2 (3) 

[(27r) M + 1 detC(^)] 1/2 

has been shown in [T7] to describe well the statistics of eigenstates of perturbed cat maps 
and to lead to a quantitative understanding of the influence of scarring on the statistics of 
tunnelling from chaotic potential wells. The distribution has been written for a finite number 
of components (xq, x±, • ■ ■ , Xm)- The statistics of probe states sufficiently localised near the 
periodic orbit (but otherwise arbitrary) can straightforwardly be obtained from the limit M — > 
oo. More generally though, this limit needs to be treated carefully, as discussed in Section 2. 

In this paper we investigate the local structure of scarred eigenfunctions by using the enve- 
lope matrix C(9) to define a series of probe states adapted to a given periodic orbit and a given 
part of the spectrum. These scar states are defined so that the joint-probability distribution 
P(x; 9) predicts statistical independence for their overlaps with chaotic wavef unctions. In a 
phase-space representation, these probe states are concentrated along the stable and unstable 
manifolds of the periodic orbit. In a region of the spectrum associated with maximal scarring, 
the first of these states corresponds to the scarfunction introduced by Vergini and Carlo in 
[TH - Elsewhere they correspond to linearised eigenstates such as considered by Nonnenmacher 
and Voros ^3] and Kaplan and Heller jHj. The analysis here suggests an extension to a series 
of such states, however, forming a local basis near the periodic orbit. Furthermore these scar 
states vary as a function of the spectral parameter so that the local basis varies quasiperiodi- 
cally as the spectrum is traversed. A surprising outcome, which has previously been noted in 
[U], is that in antiscarred regions of the spectrum, where conventionally we expect the chaotic 
wavefunction to be suppressed in a neighbourhood of the periodic orbit, we can define probe 
states for which chaotic wavefunctions have larger overlaps than predicted by RMT. That is, 
we find positive scarring in the antiscarred regions. We see in this paper that this observation 
generalises to a basis of such states, in terms of which a statistical description of wavefunction 
overlaps can be achieved using a simple product of independent Gaussians. 

The organisation of the paper is as follows. In Section 2 we start by recalling the definition 
of the envelope matrix C{9) and then we define optimised probe states which give decoupled 
overlap distributions with scarred eigenfunctions. We also discuss the local pattern of the 
scarred eigenfunctions and how these vary with the spectral parameter 9. We see how this 
works for a specific quantum map model in Section 3 and we consider billiard systems in Section 
4, where we show that the statistics of the boundary eigenfunctions can also be described by 
the joint-probability distribution P(x; 9). 



2 The envelope matrix and scar states 

In this section we begin with a brief description of the envelope matrix C(9) which determines 
the joint-probability distribution P(x; 9). The eigenvectors of C(9) can be interpreted as pro- 
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viding the statistically independent basis for the joint-probability distribution. The eigenvectors 
with the largest eigenvalues provide us with probe states which give overlap distributions in 
which the effect of scarring is optimised. Borrowing the terminology of Vergini and Carlo, we 
refer to these probe states as scar states. We show in particular that the first of these scar states 
coincides in the semiclassical limit with the scarfunction introduced by Vergini and Carlo ^1] 
and that in general they form a subset of the hyperbolic stationary states whose coherent-state 
representations have been described in detail by Nonnenmacher and Voros [T3] . 

The notation in this section is adapted to the case of a quantum map U with eigensolutions 
U\ip n ) = e~ ien \tp n ). We consider statistics of ensembles which are constructed by restricting 
the eigenangle 9 n to a narrow window centred at the value 9. The corresponding construction 
in the case of time-independent systems is to consider states for which the action S(E n ) of the 
periodic orbit under consideration is such that 6 n = S(E n )/h is similarly restricted. 

2.1 Envelope matrix C(#) 

The joint-probability distribution P(x; 8) completely characterises the influence of scarring 
on wavefunction statistics near the periodic orbit and our assertion is that this is in turn 
completely determined (to a good approximation at least) by the envelope matrix C(8). For 
a given periodic orbit, the matrix C(6) is constructed in practice from a Fourier transform 
of a correlation matrix A(t), formed by writing linearised quantum evolution in a harmonic- 
oscillator basis For simplicity of notation, let us assume that the periodic orbit is in 
fact a fixed point p of a chaotic map and that this map is quantised by the unitary operator U, 
while Un n quantises the linearised dynamics around Pq. The matrix A(t) is simply U* in written 
in the basis with elements 

A kl (t) = ( Vk \UL\<Pi)- (4) 

We might of course use any basis of states localised in phase space around po but we use a 
harmonic oscillator basis because we can then write closed-form expressions for Aki(t). These 
are given explicitly in Appendix El but here we simply note the essential features. 

The linearised evolution underlying U\ in has inversion through p as a symmetry (even if 
the full map U does not) and this leads to a decoupling of A(t) into blocks corresponding to 
even and odd values of k. That is, Aki(t) = unless k and I are both even or both odd. In any 
discussion of statistics it is therefore natural to consider separately the components for even 
and odd k — any explicit calculation we give will be for the even case. The linearised evolution 
also has a time- reversal symmetry (again, even if the full map U does not). This means that 
we can choose phases of \<pk) so that A(t) is symmetric and (A(t))* = A(— i). 

For a map with one degree of freedom, two classical parameters which characterise the 
classical situation are sufficient to determine A(t). The first of these is the stability eigenvalue 
A* = ±e pt of the unstable dynamics around po. This sets the rate at which A(t) decays with t 
and subsequently the strength of deviations due to scarring. The second parameter effectively 
characterises the harmonic-oscillator basis. Any elliptic linear evolution which has po as a 
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fixed point can be used to generate a harmonic-oscillator basis \<pk), and we should expect 
A(i) to depend on the eccentricity and orientation of this elliptic evolution (relative to the 
hyperbolic structure of the unstable dynamics around p ). It turns out that a single parameter, 
which we denote Q following Kaplan and Heller, is sufficient to completely characterise this 
dependence, and is easily interpreted geometrically. Let (q,p) be canonical coordinates for 
which the harmonic oscillator defining \ipk) is a multiple of q 2 + p 2 . Then Q = cot a, where 
a is the angle in coordinates (q,p) between the stable and unstable manifolds of the chaotic 
dynamics about p . 

The expressions for A(t) given in Appendix El are simplest for the case Q = 0, which corre- 
sponds to orthogonal stable and unstable manifolds (relative to the coordinates (q,p) suggested 
by the normal form for the harmonic oscillator). If we are interested in the wavefunction statis- 
tics in isolation, then the basis | tpf.) may be chosen at our convenience and it seems natural in 
that case to ensure that Q = 0. In applications to tunnelling statistics [Ej, however, the basis 
\(Pk) is imposed by a secondary calculation and it is therefore useful to be able to treat the case 
Q 7^ 0. We will also find this option convenient when we apply the results to an explicit model 
in the next section. 

Given A(t), we define the envelope matrix by 

oo 

C(0) = ]T e m A{t). 

t=—oo 

Given the properties of A(i) described above, this is a real symmetric matrix which decouples 
into even and odd blocks. If we construct an ensemble of chaotic eigenstates \ip n ) from a small 
window around the spectral parameter 9, modulo 2tt, then it can be shown that the components 
Cik{9) give the averages (x*Xk) ^Zj. Given such an RMT- violating constraint on variance and 
correlation, a common procedure in the analysis of wavefunction statistics has been to assert 
that the statistical distribution remains normal but with an appropriately adjusted covariance 
matrix. In particular, this has been done by Antonsen et al [TH] for averages of the wavefunction 
along the periodic orbit and by Narimanov et al [3] for the wavefunction itself. In our case this 
amounts to assuming the joint probability distribution (jBJ) in the case of GOE statistics, with 
an obvious generalisation for the GUE case. While we cannot prove the result, the distribution 
was shown in to successfully describe the statistics of low-lying components of x and to 
explain deviations due to scarring of tunnelling-rate statistics. We posit the result simply as a 
conjecture and refer to the numerical evidence presented here and in jT2j as justification for it. 

The joint-probability distribution obtained from C(6) is written in (J3J) for a truncated set 
of components (x = xq, xi, ■ ■ ■ , Xm)- To obtain a complete characterisation of wavefunction 
statistics, we should take a limit M — > oo. However, because C(6) is derived from a linearised 
hyperbolic system whose quantisation does not have normalised eigenstates, the elements of 
C(8) decay quite slowly and the limit must be treated with care. We find a straightforward 
limit if we use (jHJ) to obtain the statistics of a measure such as 

U = (V>„|P|Vn) (5) 
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which samples the chaotic state \ip n ) in a localised region of phase space around po- Localisa- 
tion in the context of © means that P is an operator whose matrix representation P in the 
harmonic-oscillator basis decays sufficiently rapidly that TrPC(#) exists — it has been shown 
in J7j that the statistics of £ n are then governed by a distribution of the form (assuming GOE 
statistics and Hermitian P) 

P(6 *) = 7T / ; , = dg, (6) 

2nJ v /det(l + 2igPC(0)) 



that is well-defined in the limit M — > oo. Alternatively, any finite subset of the components 
x = (xo, xi, • • •) can be described by a joint-probability distribution written analogously to © 
with an appropriate submatrix of C(6). The limit M — > oo is not simple however, if we try to 
treat P(x; #) as a full distribution. 



2.2 The construction of scar states 

Consider the distribution for a finite set of coordinates x = (x , xi,"', %m)- As for any normal 
distribution it is natural to rotate these coordinates so as to produce a statistically independent 
basis y = (y , y u ■ ■ ■ , y M ) for which 

(ViVj) = Si{e)5ij (7) 
and for which the joint-probability distribution is a simple product of Gaussian distributions, 

M M i 

p(t, o) = U P(y>, e) = Yl e-yy^ e \ (8) 

i=o i=o ^2-KSi(6) 

These are obtained as the eigenvectors and eigenvalues of the corresponding (M+ 1) x (M+ 1) 
submatrix of C(#). It should be noted, however, that because the infinite-dimensional limit of 
C(6>) does not have proper eigenvectors, the coordinates yi and the eigenvalues 

s (e) = (s (9), Sl (e),---,s M (e)), 

do not converge in a simple way as M — > oo. We will see in particular that the leading 
eigenvalues diverge logarithmically with increasing M. This means that we can construct 
linear combinations of the harmonic-oscillator basis states \(pk) for which the overlaps with the 
chaotic states \ip n ) can be made arbitrarily large, on average, compared with the expectation 
from RMT. It should be noted however that for any finite value of Ti there is a limit to how large 
we can let M become without the linear approximation underlying the construction of C(8) 
becoming invalid and the limit M — > oo can therefore only sensibly be taken in conjunction 
with the semiclassical limit [T§] . 

We begin by deriving formally the eigenvectors of C(8) in an infinite basis and later indicate 
how these solutions are regularised by truncation. For this purpose we use the fact that A(i) 
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and C(9) decompose into even and odd blocks and consider these blocks, denoted by A ± (t) 
and C ± (9), separately. Using Poisson resummation we write 

00 POO 

C^O) = / e ii9+2nm)t A ± {t)dt 

m=-oo 
00 

= £ F ± (# + 2vrm) (9) 



and we claim that 



/oo 
e int A ± (t) dt 
-00 



is an unnormalised projection onto an eigenvector of C ± (^) if Q = 9 + 2-nm for some integer 
m. We see this by writing 

F±(n) = ( Vk \S(Q)\Vi), 

where the superscript ± indicates whether \(pk) and \(fi) are restricted to the even or odd 
subspaces and S(Q) is the operator 

POO 

s(n) = / e int Ul in dt. 

J —00 

Let us now substitute 

fjt _ p -iHt/n 

U ]in — e 

where H is a quadratic Hamiltonian linearising the hyberbolic dynamics about pq. We expand 

^nn = E / dEe-^lE^iE^ 
± J 

in eigenbasis HlE^) = E\E^) of H, with the usual normalisation (E ± \E' ± ) = 5(E — E'). The 
sum over ± here arises because the eigenstates for each E are doubly-degenerate. In a repre- 
sentation where the problem amounts to scattering from a parabolic barrier, the conventional 
description of this degeneracy is in terms of scattering states with left or right incoming waves, 
but in our case we choose eigenstates that have even or odd parity. This decomposition gives 
in turn that 

S{Q) = 2n5(Q - H/h) = 27rhY,\E{ty ± )(E{tt) ± \ 

± 

where E(Q) = HQ. We note that phase space representations of the states {E^ have been 
analysed in detail by Nonnenmacher and Voros (beware however that in their notation the 
label ± sometimes refers to the direction of an incoming wave whereas here it always refers 
to parity) and we will return to their phase-space description later. This formal result states 
simply that the eigenvectors of C ± (8) are representations in the harmonic-oscillator basis of 
the subset of these states for which Q = 9 + 27cm and which have a definite parity. In particular 
we may decompose C ± (^) as promised in 0, with 

F ± {Q) = 2ixh w ± {Q) w ± (fi) t (10) 
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and w£(Sl) = (ifklEin)*). 

Equation (jHJ) therefore provides us formally with a spectral decomposition of C ± (6) with 
eigenvectors w ± (6 l + 2irm) labelled by the integers m. The corresponding eigenvalues diverge, 
however, because the hyperbolic stationary states |-E' ± ) are unnormalisable. This supports our 
earlier assertion that the eigenvalues of a finite submatrix of C(9) diverge as M — > oo. To 
regularise the situation we must consider truncations of C(6). This corresponds physically to 
restricting the probe states we use to be sufficiently localised near p Q in phase space, which is in 
any case a necessity if the linearised evolution used in deriving the form of A(t) is to be valid. 

We will consider truncations of the form 

C kl (0) = y/P~ k C kl (0)^Pi 

where the coefficients P k decay from unity sufficiently rapidly that C(9) has a series of non- 
vanishing eigenvalues with normalisable eigenvectors, but sufficiently slowly that C(9) can be 
regarded as an approximation for C(9) in an appropriate neighbourhood of pq. In the simplest 
case we use a step function Pi = 0(M — /), which replaces C(9) with a finite submatrix describ- 
ing the statistics of (xo,xi, ■ ■ ■ ,xm), but we may also let P k decay smoothly. This truncation 
is useful in the statistics of (0) above in the case P = J2k Pkl^kjifk] (which is the form seen in 
applications to tunnelling, where P k is a tunnelling rate associated with the probe state \<fk))- 
In particular we will use truncating functions of the form 

Pk = 1 ZjZ e P(k-M) ' ( n ) 

where M behaves as a cutoff in the dimension of C(9) and j3 provides a decay rate after cutoff. 
The sharp truncation above corresponds to the limit (3 — > oo. 

The eigenvectors of the truncated matrix C(6) are not as straightforwardly calculated as 
those of the full matrix C(6). The leading eigenvectors, however, corresponding to smaller 
absolute values of m, should approximately coincide with those of C(8). In detail, if E(Q) with 
Q = 9 + 27rm is small enough and the components of P decay slowly enough from unity then 

(^(0^ = ^(0) 

is approximately an eigenvector of C(9) with eigenvalue 

s(0) = 27r^(v / Pw ± (fi)) t v / Pw ± (fi) 

= 27ih(E(Q) ± \P\E{Q) ± ). 

Since ^(f])*) is unnormalisable, these eigenvalues diverge as promised as the truncation di- 
mension increases and P — > /. The eigensolution is only approximate because, while (|10|) still 
applies to the truncated matrix if w ± (f2) is replaced by V / Pw ± (fi), the vectors \/Pw ± (r2) for 
two different values of Q are only approximately orthogonal. 

We present in Figure [T] the leading eigenvalues, relabelled (so(9), si(9), ■ • •) in order of de- 
creasing value, for a fixed point po with stable and unstable manifolds along p = ±q, a harmonic 
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Figure 1: The leading eigenvalues s(9) = (s Q (9), Si(9), • • •) of the envelope matrix C(9). Here we 
use truncation parameters (3 = 0.025 and M = 300. The bullets indicate Si(0) for % = 0, 1 and 
2 and for 9 = —n/2, 0, 7r/2 and 7r. The Husimi distributions corresponding to these eigenvalues 
are shown in Figure 121 

oscillator basis with Q = and using a smooth truncation of the form given in ([lip. As ex- 
pected from the preceding discussion, these eigenvalues unfold to give a smooth function of 
Q = 9 + 27rm, decreasing from the centre. The leading eigenvalue is peaked in the so-called 
scarred region of the spectrum — we have chosen phases so that this occurs at 9 = 0. In terms 
of wavefunction statistics, this means that the corresponding eigenvector defines a probe state 
for which the variance is much larger than predicted by RMT (about 18 times larger for the 
parameters illustrated in Figure^). The scarred region is also where overlap statistics of chaotic 
wavefunctions with a simple Gaussian probe state have the greatest variance. In the case of a 
Gaussian probe state, however, the variance falls between these peaks to a value significantly 
smaller than predicted by RMT in the anti-scarred region, corresponding here to the region 
around 9 = n — this contrasts with the case where we track a probe state corresponding to 
the leading eigenvector, in which case the overlap statistics have large variances even in the 
anti-scarred region (albeit not as large as in the scarred region). 

The leading eigenvalue sq(9) can be tracked through crossings in the antiscarred region to 
give the next eigenvalue Si(9) and so on. Although less dramatically deviant than the leading 
eigenvalue, these subsequent eigenvalues define probe states which also have overlaps that are 
significantly larger than predicted by RMT. As with the leading probe state, the overlaps they 
define are larger than RMT even in the antiscarred region. For any finite truncation dimension 
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and increasing index, the eigenvalues will eventually become very sensitive to the truncation 
and will not define useful probe states. For a fixed index and increasing truncation dimension, 
however, the probes states provide us with a useful characterisation of wavefunctions near the 
periodic orbit and we will refer to these probe states as scar states. 

We will see later that the leading eigenvalues increase logarithmically with the effective 
truncation dimension M. The leading eigenvectors, on the other hand, converge to give a 
meaningful limit in the region of phase space surrounding the periodic orbit. A detailed dis- 
cussion for the Husimi distributions of such hyperbolic stationary states has been given in [T3] 
and we illustrate these with some explicit cases here. Given an eigenvalue Si(9), we denote the 
corresponding scar state by \Xi{9)) which in the limit M — > oo approaches 



Figure 121 illustrates this Husimi distribution for some of the leading scar states, whose eigen- 
values Si(9) are indicated by bullets in Figured The scar states shown are appropriate to 
a fixed point po at the origin of phase space with stable and unstable manifolds defined by 
p = ±q. Each of these Husimi functions is strongly localised around the stable and unstable 
manifolds of the fixed point, although the details vary with the spectral parameter 9 and index 
i. In distinguishing between different cases it is useful to note that the zeros of the Husimi 
functions are visible as dark spots — these enable us to distinguish between states which are 
concentrated in similar regions of phase space but which have different phase structure. 

The leading scar states \Xi(9)) are potentially useful not only for understanding the statistics 
of wavefunctions but also for characterising the structure of individual wavefunctions. The 
variance of the overlap (xiifylrfn) is quite large and in fact can diverge in the semiclassical 
limit. In a given part of the spectrum, chaotic eigenfunctions will therefore tend to have 
especially large components along the corresponding scar states. One aspect of this is that 
we might therefore recognise the characteristics of the scar states in individual states, and the 
changing form of the scar states across the spectrum will provide a template for recognising such 
structures. We will see explicit examples of this in the next section. A second aspect is that 
the scar states should provide a powerful guide to choosing bases for the explicit calculation of 
chaotic eigenstates. This sort of construction has recently been proposed by Vergini and Carlo 
PU] . They have in particular defined for each periodic orbit a scarfunction which is concentrated 
around the stable and unstable manifolds much like the scar states \Xi{9)) are. In fact we will 
see in the next subsection that the scarfunction coincides with the leading state |xo(0)) in the 
semiclassical limit. The property of the states \xo(9)) that they change as a function of 9 and 
adapt to a given part of the spectrum, as well as the ability to generate additional orthogonal 
states (labelled by i), might therefore provide a useful generalisation of that approach. 



Xi {9))^\E{9 + 2im H ) ± ) 



for some integer va^ We denote the corresponding Husimi distribution by 
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Figure 2: The Husimi functions of the scar states corresponding to the bullets in Figure Q The 
underlying hyperbolic system has stable and unstable manifolds along the diagonals p = ±q. 
The first row gives Ho(q,p,9) for 9 = — 7r/2,0,7r/2 and 7r. The second and third rows give 
Hi(q,p,9) and H2(q,p,0) for the same values of 9. Zeros of the Husimi functions appear as 
dark points. 

2.3 The leading scar state 

We show now how the decomposition outlined in the previous subsection emerges in the 
harmonic-oscillator basis. We cannot give closed-form expressions for general scar states \Xi(9)) 
in this basis, but we can calculate the basic properties of the leading state |xo(0)) in the scarred 
part of the spectrum, corresponding to 9 = in our convention. We will see in particular that 
this reproduces the scarfunction of Vergini and Carlo . 

We have noted that the envelope matrix C(9) depends on the stability exponent p of the 
underlying hyperbolic dynamics and a parameter Q depending on the basis. We are free to 
choose a basis for which Q = and we assume that to be the case in this subsection. The 
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stability exponent p determines how quickly the correlation matrix A(t) decays with t and 
therefore how far the envelope matrix C(8) deviates from the RMT case C(8) = I. The 
eigenvalues s(6) in particular depend strongly on p. The scar states, however, depend only on 
the underlying geometry of the stable and unstable manifolds and not on the rate of decay. We 
therefore consider the limit p — > 0, which simplifies because we can then replace the Fourier 
series defining C(9) by an integral 

J2Mt)e m ^ [ A(t)e i9t dt. (13) 
t •* 

In other words, the sum in (JUJ) is dominated by the m = term. We note that this formulation 
is similar to the construction of the universal scar measure in jH] or the continuous quasimode in 
[TE] , except here we also incorporate the propagation of excited harmonic-oscillator eigenstates 
as well as the Gaussian ground state. 

Using the form for A(t) given in Appendix [X] we can in particular evaluate this integral 
explicitly for the case 6 = 0. We find a decomposition of the form 

C(0) ~ -vv T , (14) 
P 

where the components of the vector v are 



, ir(27rrl /4 y'N r(n + i) = 

j 1 ] 22«(2n)!i> + !) ° K n ' 



(15) 
otherwise. 



This means that v is the leading eigenvector of C(0) and that the corresponding eigenvalue is 

s (0) ~ iv T v. (16) 
P 

The vanishing of components for odd k is a reflection of the fact that linearised evolution has 
inversion through the origin as a symmetry — we have already see that this symmetry leads 
in the general case to a decoupling of C(9) into odd and even blocks. When we choose a basis 
for which Q = 0, the stable and unstable manifolds are orthogonal and we have in addition 
a symmetry of interchange (q, p) — > (p, —q) of the corresponding canonical coordinates, this 
is responsible for a further decoupling of the index k mod 4. This four-fold symmetry is also 
reflected in the general case, in particular through the symmetry Hi(q,p,9) = Hi(p,q, —6) 
visible in the Husimi functions of Figure El 

We remark that this eigensolution depends on the stability exponent p only through the 
prefactor in pfij) . The calculation has been performed formally for the untruncated matrix 
C(8) and we note that (fTH|) is divergent and gives a meaningful eigenvalue only if we truncate 
the matrix at a finite dimension M. We can obtain an estimate for the divergence of Sq(0) with 
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truncation dimension by noting that the components of v have the asymptotic behaviour (as 
k = An — > oo) 

2 1 

^4n - n 

The leading eigenvalue for a matrix sharply truncated at dimension M therefore asymptotes to 
(as p — > and M — ■> oo) 

1 M/4 1 

s (0) ~ - V vL ~ - log(M/4) + const. (17) 

Although our derivation has been restricted to the limit p — > 0, we remark that this describes 
quite well the divergence of sq(0) even for moderate values of p, as indicated in Figure El where 
we compare it with a numerical calculation for p = 0.6. 
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10 10 2 

M 

Figure 3: The leading eigenvalue Sq(0) as a function of the sharp cut-off dimension M when 
p = 0.6 and Q = 0, showing a logarithmic divergence. The solid curve is the approximation in 

(ED 

We expect similar divergence to hold for the next-to- leading eigenvalues. Since Sj(0) di- 
verges we might then expect to be able to make the statistical deviation from RMT arbitrarily 
large by choosing appropriate probe states. It should be noted, however, that our results are 
meaningful only if we take the limit M — > oo in conjunction with the semiclassical limit H —>■ 0. 
Our calculation of the correlation Cui{0) between two basis states was obtained by assuming 
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evolution that is linearised about the periodic orbit. This is always valid for fixed k and I in the 
semiclassical limit. If we fix h and take the limit k, I — » oo, however, the basis states \(pk) will 
eventually fall outside the region of phase space where the linearisation is appropriate and the 
correlations Cki{0) will then not accurately describe the statistics. For any given quantum sys- 
tem then, we should restrict ourselves to truncations for which linearised evolution is accurate 
and this effectively limits the value of Si(9) and therefore the deviation from RMT. 

Consistent with the preceding subsection, the vector v can also be expressed formally as an 
eigenstate of the Hamiltonian 



H 



P 



(P 2 ~ q 2 



giving the underlying hyperbolic motion in appropriate canonical coordinates (q,p). We can 
express this Hamiltonian in terms of creation and annihilation operators of a harmonic oscillator 
basis for which Q = as 



H 



P 



a' + a 



The vector v is then easily seen formally to be a null eigenvector 



Hv = 



(19) 



of the Hamiltonian matrix 



(H) w = (<p k \H\<pt). 

In particular the leading scar state is also one of the simpler cases for which Husimi functions 
were constructed in [13] . This can be seen by using the representation 

1 



-1*1=72 



k=0 



\fh) 



of the coherent states, where z = (q — ip)/ \/2h, and substituting the coefficients of ([T5]) in (p"2 
giving 

2 



H o (z,0) = -Le- ML ' 



n=0 1 V L ' 4 



v 4n 



r(n+ f)2 2 «(2n)! 



2V 2 



z 



\fzJ_iu — 



for the Husimi function of this state. This coincides with the expression in [J3] for the E = 
eigenstate of the hyperbolic system as expected from the discussion of the previous subsection. 

The preceding discussion also helps us to relate |xo(0)) to the scarfunction of Vergini and 
Carlo [Tl] . They define the scar function \x) to be a finite linear combination of the states \<pk) 
which minimises 

(X\H 2 \ X ). 

In the limit M — > the eigenvector v provides a solution to just this problem and the corre- 
sponding state |xo(0)) should therefore coincide with the scarfunction of [Jl] in that limit. 
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Away from the scarred region, that is for 9^0, the scar states should correspond to universal 
test states [El HE], or quasimodes [IB], although we have not written explicit forms for them in a 
harmonic-oscillator representation. Note however that one can find explicit expression for the 
Husimi functions in ^3]. We also note that phase space representations of operators with the 
structure of C(6) have also been given by Rivas and Ozorio de Almeida [21], and these might 
also be useful in deriving explicit expressions for the scar states in phase space. 

3 Application to a quantum map model 

We now investigate how well the scar states describe wavefunction structure and statistics for a 
specific quantum map model. The details of this model are unimportant and we merely remark 
that the classical system is a hyperbolic mapping of the torus with time-reversal symmetry and 
a fixed point at the origin — the detailed parameters for the system are the same as those 
given in and the reader should look there for further details. As in |17j . we choose as 
basis states \<pk) the eigenstates of a Harper Hamiltonian which is harmonic near the fixed 
point. This problem is characterised by a stability exponent p = 0.6006 and a basis-parameter 
Q = 0.5658. The semiclassical limit is controlled by the dimension of Hilbert space, which we 
denote by N. Because of numerical limitations we use relatively modest truncations of C(9) in 
this investigation — typically we use a sharp truncation with M = 28 or less. 

We begin by considering some individual wavef unctions. We show in FigurelUsome examples 
of eigenfunctions of the quantum perturbed cat map which reflect the structure of the leading 
scar state in Figure El at the corresponding eigenangle. The enhancement of the eigenstates 
along the stable and unstable manifolds as seen here is a well-known characteristic of scarred 
eigenfunctions. However, the dependence of the detailed pattern of this enhancement on eige- 
nangle and the fact that it can persist into the antiscarred region, is less obvious. Although we 
do not show examples here, we also expect to find eigenfunctions with enhanced overlaps with 
subsequent scar states with % = 1, 2, • • •, whose Husimi functions vary across the spectrum as 
indicated in Figure El The Schnirelman theorem [22] guarantees that such structures cannot 
dominate individual wavefunctions in the semiclassical limit, However we note that the scar 
states shrink in phase space as h —>■ and may dominate in a sufficiently small neighbourhood. 
We also note that even when scar states do not dominate a given wavefunction, locally or oth- 
erwise, they should collectively form an efficient basis for calculating and characterising such 
wavefunctions [2T)] . 

Let us now turn to overlap statistics. We have seen that if we use the scar states \Xi(9)) as 
probe states, the statistics of the resulting overlaps 

l/i=(tt(Wn>, (20) 

should be described by the Gaussian distribution P(yi] 9) given in Eq. (JBJ). Furthermore, over- 
laps for different values of i should be statistically independent. We now test this supposition 
using the perturbed cat-map model. 
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Figure 4: The Husimi functions of some scarred eigenstates of the perturbed cat map. Starting 
from the left, these have the Hilbert-space dimensions and eigenangles (N = 404, 9 = —1.690), 
(N = 400,0 = 0.050), (N = 400,0 = 1.547) and (N = 404,0 = 2.995). We note that 
eigenangles are respectively close to — 7r/2,0,7r/2 and 7r and that the states bear a qualitative 
resemblance to the leading scar function at corresponding values of 9. 

To verify statistical independence of distinct components y^ we consider the variable 

£ = Vi Vj with % ^ j. 
The distribution for this variable can be calculated as a special case of (JBJ), to give 



TzJ Si {9) Sj {9) 




Si{0)s J {9) i 



(21) 



where K (z) is a modified Bessel function. The independence of yi and jjj is reflected in the 
fact that this is an even function of £. Figure verifies that distribution describes the joint 
statistics of yo and y± in the case of a perturbed cat map and supports the assertion that these 
variables are statistically independent. Note that numerical limitations restrict us to relatively 
modest truncations of the harmonic-oscillator basis in this case — the statistics shown have 
been computed for scar states with M = 28 and (3 = oo — nevertheless the agreement is quite 
good. 
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Figure 5: The distribution of £ = yoyi, where y- t = (Xi(O)lV'n)- The histogram shows the 
results of a numerical calculation for a perturbed cat map in which the eigenstates \ip n ) in the 
interval — 0.057T < 9 < 0.057T used. The solid line is the theoretical expectation Eq. (|21jl . The 
symmetry of this result supports the assertion of statistical independence for y Q and y±. 

We also consider the statistics of a single component, concentrating on the statistics of 
yo, the overlap with the leading scar state. We first construct the probe state |xo(0)) for the 
maximally scarred part of the spectrum and examine its overlap statistics as a function of 9. 
The variance at any other part of the spectrum is then calculated for this fixed state using 

CJfi) = E^(Xo(0)|^ n |Xo(0)). (22) 

t 

This is analogous to the envelope for the universal test state considered in j^J. A comparison 
with numerical results for the perturbed cat map is shown in Figure El for two truncation 
dimensions M of the matrix C(9). We see that the overlaps are concentrated more in the 
scarred region and that the height increases as M is increased, consistent with the observation 
that s (0) diverges logarithmically with M. Note that, as for any fixed state, a spectral average 
of the variance must average to unity since (C sc ) = (xo(Q)\U^°\xo{^)) = 1. Any growth in the 
scarred region must therefore be accompanied by a decay elsewhere in the spectrum and this 
is also visible in Figure |H1 This does not contradict the fact that Sq(6) remains large over the 
whole spectrum because the leading scar state varies with 9 whereas in this calculation we have 
used a fixed state |xq(0)) 

We also show in Figure Efc) and (d) the distribution of y — (xo(O)lV'n) about this mean 
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Figure 6: The variance calculated for the scar state |xo(0)) according to ()22)1 is compared with 
results for the perturbed cat map as a function of 9. Cases (a) and (b) are respectively for 
effective truncation dimensions M = 8 and M = 28 with (3 = oo in the scarred region. The 
accumulated overlap distributions P(y) are shown in (c) and (d), respectively, for the same 
truncation dimensions. The solid line comes from the present theory and the dotted line is the 
RMT prediction. 



behaviour, averaged over eigenangle 6. The distribution 



P(y) = ^ / d8 P (y, 9) = i / d9 1 e »V**.W (23 ) 



shown as a solid line describes quite accurately the wavefunction statistics, shown as a his- 
togram. Note the marked difference between these results and the RMT prediction indicated 
by a dotted line. Similar reults are found for the statistics of (xoifyl^n) with i > 0, although 
in the deviation from RMT is then less stark. 
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4 Wavefunction statistics in billiard systems 



The discussion so far has been in terms of quantum maps, but we now address the physically 
more interesting question of time-independent systems, and in particular billiard systems. In 
this section we describe how the the joint-probability distribution may be applied to such 
systems, for which an analogy with maps can be made through the boundary-integral formalism. 
In more generic problems such a connection is made using the transfer-operator formalism of 
Bogomolny [23] . 



4.1 Boundary eigenfunctions 

We consider eigenfunctions of the Helmholtz equation 

(V 2 + A£)* n (x) = 0, 



(24) 



satisfying the Dirichlet condition \I/„(x) = on the boundary and we assume that the domain 
is such that the corresponding billiard dynamics is chaotic. In order to apply the results of the 
preceding sections it is convenient to reformulate this problem in terms of maps, and a natural 
means of doing this is to use the boundary-integral method, which provides a quantum analog 
of the classical Poincare-Birkhoff mapping. 

The boundary- integral method has been extensively described elsewhere so here we just 
describe the key features, primarily to establish notation. For a more detailed description, and 
in particular for other investigations of scarring in the boundary representations, we refer to 
[2U I2EJ EE1 123 EEl I2E|- The essential point is that we replace the Helmholtz equation in the 
interior of the billiard by the integral equation 



ip n (s) = / K(s,s';k n )ip n (s')ds' 



over a coordinate s on the boundary, where 



dn 



(25) 



(26) 



x=x(s) 



serves as a boundary eigenfunction. The kernel is 

<9G(x, x(s'); k n 



K(s,s';k n ) 



dn 



x=x(s) 



where 



G(x,x';A ;n ) = ^ 1) (A; n |x-x / | 



(27) 



(2f 



is the Green function for the Helmholtz equation. 

The integral equation (|2*K|) will play the role of a quantum map for us. Even though the map 
is not unitary [24J, we will assert that a straightforward generalisation of the envelope matrix 
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derived in Appendix A constrains the variances of overlaps between harmonic probe states and 
eigenfunctions. The important point is that iterates of the map concatenate semiclassically in 
the same way that iterates of a unitary map do and while there are some local modifications of 
amplitude on classical length scales due to breaking of unitarity, these do not appreciably affect 
the statistics when phase-space localised probe states are used. In particular, we assert that 
the statistics of if) n (s) near a periodic orbit are described by the joint-probability distribution 
with the structure described in previous sections, at least in the semiclassical limit. In order to 
write this distribution down we need to define analogues for some of the ingredients we used in 
the previous sections. 

In describing the Poincare-Birkhoff mapping we use canonical coordinates (s, v) where s 
is an arc-length coordinate at which a trajectory collides with the boundary and v = sin a, 
where a is the angle of incidence. We consider wavefunction statistics in the neighbourhood of 
a point with coordinates (so, Vq) where a periodic orbit of length £ intersects the boundary. As 
a spectral parameter we use 9 = k£ — /ivr/2, where /i is the Maslov index of the periodic orbit 
- with this definition we will find that that full scarring occurs when 6 = mod 2ir. As probe 
states we use eigenstates of a harmonic oscillator centred at (so,fo) on the boundary section. 
These are written as boundary wavefunctions of the form 



1 

<M S J = , exp 

f(T2 l l\y/7r 



2a 



+ ikv (s — s ) 



— (29) 



a 



where Hi(z) is a Hermite polynomial. As in [7j we let a, which provides us with an aspect ratio 
in phase space of the ellipse defined by the harmonic oscillator, to scale with k as 



a 



where £q is a characteristic length of the billiard. In practice we choose £q so that Q = (that 
is, so that the invariant manifolds are orthogonal in the corresponding metric). We scale the 
overlap variables x = (xq, Xi, ■ ■ •) as 

xi ~ \Zk{cpi\ip n ) = \fk I ds cp*(s)ijj n (s), (30) 



o 

where L is the boundary arc length. These scalings are such that the statistics of the components 
Xi do not vary over the spectrum. We also note that kL plays the same role as the dimension 
iV in the case of quantum map. 

For numerical purposes we consider the cardioid billiard, which is defined in polar coordi- 
nates by 

r = l + cos#. (31) 

We consider just the odd-parity states and therefore restrict ourselves to the half-cardioid 
illustrated in Figure0 For simplicity we treat the shortest periodic orbit, also shown in Figure[3 
which has stability exponent p ~ 0.693 and length £ = 3\/3/2. We examine wavefunction 
statistics where this intersects the boundary section at sq ~ 2.00 and vq = 0. 
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Figure 7: The cardioid billiard and the periodic orbit used in numerical calculation. 
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Figure 8: The overlap probabilities y\ = k\ (xol^n) | 2 for the leading scar state are shown as a 
function of k. In (a) the state xo is constructed for the scarred part of the spectrum with 8 = 
and in (b) it is constructed for the antiscarred part, with 6 = ir. In each case the scar state has 
been calculated using a truncation dimension of M — 12. 



We calculate scar states as described in the previous sections, with the obvious modifications 
appropriate to billiards, and denote the corresponding boundary functions by Xi( s '^)- Let us 
begin by considering the overlap probability 

yl = k\( Xo (OMn)\ 2 , (32) 
for the leading scar state calculated for the scarred part of the spectrum. The scaling with k 
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here is designed so that the average value does not change across the spectrum. This overlap 
probability is calculated for individual eigenfunctions and the results shown in Figure IBfa) as 
a function of k. As expected we find strong peaks in typical values localised around values 
of k for which 9 = ki — fnr/2 = mod 2tt, with a compensating suppression between peaks. 
The sharp contrast between the peaks and valleys in this figure is a reflection of the extent to 
which xo(<s;0) maximises the deviation of statistics from RMT. We also show, in Figure El^b) 
the overlap probabilities for the leading scar state Xo( s ; 7r) constructed for the antiscarred part 
of the spectrum where 9 = tt mod 2tt. Once again we find periodic peaks in the average value 
of the overlap probability, albeit with a somewhat reduced contrast. Note also, however, that 
the peaks in Figure Efb) occur when there are troughs in Figure Efa) and vice versa. In other 
words, we have defined a probe state for which significantly larger than average overlaps occur 
in parts of the spectrum which are normally associated with antiscarring. Although we do not 
show them explicitly, each of these sequences can be defined also for subsequent scar states 
Xi{s; 9) with i > 0; these subsequent sequences peak in the same parts of the spectrum but are 
statistically independent. 

To make a more quantitative comparison with the theory of the scar state we collate the 
overlap probabilities as a function of 9 and compare with the prediction of (|2*2*|) (and its general- 
isation to the antiscarred region) in Figures Ela) and (b). These describe envelopes for Figure |H1 
and confirm that the detailed predictions for the variances using the envelop matrix C(9) work 
quantitatively in this billiard system. A comparison of the statistics of the distribution of 
individual states about this average in (c) and (d) also supports the hypothesis of Gaussian 
statistics in this context. Calculations of correlation and the statistics of higher scar states 
work similarly to the case of quantum maps, as investigated in the previous section, and we do 
not give further numerical details in the boundary representation here. Instead we turn to a 
discussion of how we may formulate these statistical properties in terms of full wavefunctions 
in the interior of the billiard. 

4.2 Full eigenfunctions 

The scar states for billiards have so far been given in the boundary-integral formalism. We 
now propose corresponding probe states which are appropriate to full wavefunctions in the 
interior of the billiard. Given an eigenfunction ip n {s) of the boundary- integral equation, the 
full eigenfunction ^ n (x) is obtained from it using Green's identity, giving 



We define probe states in the interior of the billiard in analogy with this equation. That is, 
given a probe state ip(s) in the boundary representation, whether it be a harmonic-oscillator 
eigenstate, a scar state or otherwise, we define a corresponding probe state in the interior using 




(33) 




(34) 
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Figure 9: The partial variance (y 2 ) is shown as a function of 9. In (a) and (b) we respectively 
show results for the scar state of the scarred region (y 2 = k\(xo(0)\ip n )\ 2 ) and of the antiscarred 
region (y 2 = fc|(xo( 7r )|V ; n)| 2 )- The solid lines show the envelope functions C SCj o(6!) and C sc ^(9) 
calculated in analogy with (|2^j). In (c) and (d) we show the corresponding accumulated overlap 
distributions P(y) describing the distribution of individual overlaps about the averages depicted 
in (a) and (b). The solid line comes from the present theory and the dotted line is the RMT 
prediction. That is, in (c) we compare a histogram of y = k\(xo{0)\ip n )\ with theory and in (d) 
we do the same for y = k\(xo(7i)\ip n )\- 



This function automatically satisfies the Helmholtz equation in the interior but, because ip(s) 
does not satisfy the boundary integral equation, does not satisfy the Dirichlet boundary condi- 
tions placed on \l/ n (x). Nevertheless it provides a perfectly well defined function in the interior 
of the billiard and we may use it to characterise the wavefunction near the periodic orbit. 

We provide an argument in Appendix El which indicates that, given a probe state which is 
localised in phase space around a point (s ,v ), then the boundary overlap (<f\ip n ) and the full 
overlap ($|\l/ n ) approximately coincide, 

($|* ft ) ~ g{s ,po)(<p\if> n ), (35) 
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Figure 10: Examples of leading scarfunctions $Q C£ir (x; k) calculated for (clockwise from top left) 
k = 999.40, 998.80, 1000.01 and 1000.61. These correspond respectively to 9 = 0, vr/2, vr and 
3tt/2. 

up to a geometrical factor g(so,po) which depends on the location the probe state in phase 
space but is independent of ip n . Following appropriate normalisation, the predictions for the 
statistics of the boundary overlaps may therefore be carried over to the overlaps defined in the 
interior. 

In particular, we may conclude that in a given part of the spectrum the interior functions 
constructed from the leading boundary scar states Xi( s 'i ®), i — 0, 1, 2, • ■ -, should have especially 
large overlaps with the billiard eigenfunctions. Let us denote by $| car (x; k) the interior function 
obtained from a boundary scar state We show examples of these interior functions in 

Figure ITU1 corresponding to values 9 = 0, it/2, it and 3n/2 of the spectral parameter 9 and k 
near 1000. We remark that these states are somewhat extended spatially around the periodic 
orbit in configuration space, a reflection of the fact that phase space representations of the 
scar states extend quite far along the stable and unstable manifolds. We also remark that, 
near intersections of the periodic orbit with the boundary, these scarfunctions approximately 
satisfy the Dirichlet boundary conditions placed on the full solution. This is a reflection of the 
fact that they are derived from boundary states which are approximate eigensolutions of the 
boundary integral operator sufficiently near the periodic orbit in phase space. 

These functions provide us with templates for the behaviour of heavily scarred eigenfunctions 
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Figure 11: At left we show scarred eigenfunctions with, from top to bottom, k n = 92.54, 
107.48 and 118.66. These correspond respectively to 9 near (in the scarred region), —tt/2 
and 7r/2. These may be compared with the corresponding scarfunctions $Q CCtr (x; k n ), shown 
on the right. The overlap probabilities are respectively 38.2, 21.0 and 21.5 times larger than 
the RMT average. Details of the scar state can readily be detected in the eigenfunction in 
the scarred region. Although the correspondence is less obvious for the other eigenfunctions, 
features such as numbers and positions of nodes of the scarfunction can sometimes be seen in 
the eigenfunctions. 



25 



near periodic orbits in different parts of the spectrum. We compare in Figure ^2 some cases 
of individual scarred eigenf unctions with the corresponding leading scarfunction $Q Car (x; k n ). 
In some cases, such as the eigenfunction with k n = 92.54, one can identify characteristics of 
the corresponding scarfunction (such as the number of nodes seen along a periodic orbit). In 
other scarred eigenfunctions, especially those in the antiscarred region where the corresponding 
scarfunctions are not well-localised, a detailed correspondence is less evident visually, although 
there is nevertheless a large overlap with the appropriate $Q Car (x; k n ). We emphasise that, 
although not shown explicitly here, for each k n these interior functions are the leading members 
of a sequence $| car (x; k n ), % = 0, 1, 2, ■ • • and each of these functions will tend to have especially 
large overlaps with ip n . 

By tracking the scar states as a function of 9 = ki — fiTi/2 we thus define continuous 
families of interior functions which approximately vanish on the boundary near bounces of 
the periodic orbit. While our notation suggests boundary functions Xi{ s 'i@) which are 2ir- 
periodic in 9, we should note that as we increase 9 = k£ — fm/2 by 2tt, the effective value 
of h decreases and there is a slow modulation of Xi{ s 'i@) with k. (While the coefficients of 
X;(s; 0) in a basis of the harmonic states (fi(s) defined in (J29j) are asymptotically periodic in 9 
as M — > oo, the states (fii(s) themselves shrink as k increases). To emphasise this evolution, we 
label the interior functions <3>f car (x; k) with k rather than 9. As k increases by 2k /£ the form of 
$| car (x; k) transverse to the periodic orbit is approximately periodic but an extra node appears 
in the longitudinal direction. Over many such quasiperiods, there is a gradual shrinkage in the 
transverse length scale and an increase in the number of nodes along the orbit. Throughout 
this evolution, $* car (x; k) exactly satisfies the Helmholtz equation and approximately satisfies 
the Dirichlet boundary conditions near the periodic orbit. 

We remark finally that although we do not show so explicitly here, these considerations 
should also allow us to characterise scarring in general time-dependent systems. Starting with 
a formal representation of eigenstates on a Poincare section as in the Bogomolny transfer 
method, we can define an analogous extension of probe states from the section to the full 
space and therefore provide statistical predictions for overlap statistics which can be calculated 
concretely from the full wave function. 

5 Conclusion 

We have shown that the joint-probability distribution P(x; 9) describing wavefunction statistics 
in a harmonic-oscillator basis near a periodic orbit leads to the definition of a set of probe states 
which describe in a natural way the morphology of scarring in a given part of the spectrum. 
That is, in a given part of the spectrum (parametrised by 9), we define a set of probe states 
\Xii@)) with i — 0, 1, • ■ • such that the overlaps (Xi(^)|V ; n) are expected to be especially large, 
on average, and statistically independent. In particular these probe states offer a template for 
the structure of scarred wavefunctions and tell us how the shape of scarring changes across the 
spectrum. 

The scar states, in offering a basis set for which we expect enhanced overlap with chaotic 
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eigenfunctions, may in particular prove useful in explicit computation as outlined in ^3]. We 
note also the fact that we can obtain such enhanced overlap probabilities in parts of the spec- 
trum normally associated with antiscarring, where we expect a suppression of overlap with a 
Gaussian basis state on the periodic orbit. We have illustrated our results with explicit calcu- 
lation for a quantum map model and have also outlined briefly how these ideas may be applied 
to billiard systems. 
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A Appendix: The matrix elements of the linear evolu- 
tion operator 

We here give explicit formulas for the calculation of the the linear correlation matrix A(t) in a 
harmonic oscillator basis. A detailed derivation of these results has been given in J7| and here 
we state the key features. 

To begin we characterise linear classical evolution near a periodic orbit using the stability 
exponent p and the parameter Q described in the main text. We then define angles ip(t) and 
(pit) using 

coshpt + iQsinhpt = sec ijj(t) e iLp{t \ (36) 
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and we choose these angles to lie in the ranges < ip(t) < n/2 and —ip(t) < <p(t) < ip(t) 
respectively. 

Then we may calculate the matrix elements Ai k (t) for I = k + 2n and n > using the polar 
form 

A k+2n , k (t) = G kn (^(t)) e i ^ +1 / 2 )(v(*)-^)+^(t)+m,r/2 ) (37) 
where /i is the Maslov index of the periodic orbit concerned and the amplitude is 

where C%(x) denotes a Gegenbauer polynomial. We assume n > and t > in (|3T|) and use 
A* fe (— i) = A/fc(t) = to calculate the correlation function when t < or I < k. We note 

finally that A k i(t) — if k and / do not have the same parity. 

B Appendix: interior and boundary overlaps 

In this appendix we show that if a boundary probe function ip(s) is appropriately localised near 
a point (sq, vq) in phase space then there is a relationship of the form given in (|33|) between 
the boundary overlap ((p\ip n ) the full overlap ($|\I/„,) of the corresponding interior functions, 
defined by ()33|) and (j34j) . Substituting these defining relations into the interior integral for 
we find 

($|tf ft ) = f L ds f L ds' f(s,s')<p*{s)il> n tf), (38) 
Jo Jo 

where 

f(s,s') = ]f G*(x,x( s );fc)G(x,x(s');fc)dx (39) 
4 Jd 

and D denotes the billiard domain. We consider in particular the case where the probe state is 
of the form ip(s) = e thv ° s h(s) where we suppose that h(s) is supported in a neighbourhood of 
So that is small in comparison with the length scales of the billiard but varies on a scale longer 
than the length scale 1/k typical of wavef unctions. This is the case for example for the basis 
probe states ^>i(s), which vary and decay on a length scale of order \j\fk. It is also the case for 
any finite combination of them such as the truncated scar states we use. We claim that f(s, s') 
then behaves essentially as a delta-function when integrated against <p(s), giving 

J <f*(s)f(s,s , )ds~g(s , ,v )ip*(s') 

where g(s', vq) is determined from the geometry of the billiard and varies with s' on a classical 
length scale. Since (p(s) is localised to a neighbourhood of a point s that is smaller than 
classical length scales then we may further replace g(s',v ) by g(s ,Vo) in the overlap integral 
and approximate 

~ </(s ,UoX¥#n)- 
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The overlaps ($|\P n ) and ((p\ip n ) then coincide up to a factor which depends on the position 
in phase space of the probe state but is independent of ip n . This means in particular that 
predictions for the statistics of ((p\ip n ) carry over to (<&|\& n ) after appropriate normalisation. 
We justify this assertion by using the asymptotic representation of the Hankel function 



ftf>(*|x-x( a )|) ~ J— e ^|x-x( S )|- i7r /4 

V TXK x — X s 



for the Green functions in the integral (|39jh giving 

i r p t*(|x-x(s')|-|x-x(a)|) 

f(s, s')~— , dx. (40) 

8^ J d v /|x-x(s)||x-x(s / )| 

Fixing s and s' temporarily, we are free to choose the origin of coordinates to be the midpoint 
between x(s) and x(s'), so that x(s) = a/2 and x(s') = —a/2. For large the rapid oscillation in 
the integrand when s and s' are distinct means that f(s, s') is peaked around s ~ s'. To evaluate 
the integral near this peak, we consider contributions from points for which r=|x|^>a = s — s'. 
We then approximate the phase of the integrand by 

(ka i \ 

fc(|x-x(s')| - |x-x(s)|) = kasinO + O f — I , 

where (r,8) are polar coordinates for x relative the midpoint of x(s) and x(s'), which we 
approximate as lying on the boundary, and the normal to the boundary there. This substitution 
for the phase is valid as long as a <C (i? 2 /^) 1 / 3 , where R is the length scale of the billiard. We 
then approximate the integral in (|40j) by 

j « gifcasinf? 

f(s,s') ~ — - / rdrd# 

87r/c Jd r 

1 /ir/2 



r e ika a m9 L ^ de 
J-tt/2 



8nk J-tt/2 

where L(6) is the length inside the billiard of a line passing through the origin at polar angle 9. 
For convex billiards, r = L(9) simply defines the shape of the boundary in polar coordinates. 

As a function of (s, s'), this integral is effectively supported on a length scale \ jk near s ~ s' 
or a ~ 0. When integrated against a probe state which varies over a longer scale such as 1/yk, 
we may therefore treat it as a 5-like pulse. In fact, we may approximate 

V *(s)f(s,s')ds = Jh*(s)e~^ s f(s,s')ds 

•ir/2 



I roo 

-ikv 0S h *, 8 n 



8nk 



e ika( S me-v ) L ^ dQ 
w/2 



da 



l /-t/ 2 



oo 

«fca(sin6»-t) ) c | a 



e 

oo 



L(0) d^ 
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1 (p*(s') I"' 2 [2ir8(k{sine-v Q )))L{e)de 

J-tt/2 

L(e )ip*(s' 



8nk J-tt/2 
1 



Ak 2 cos 6q 



•*'-'), 



where vq = sin#o. Here, L(6q) is the length of a chord which leaves the point s' at an angle 9q 
to the normal. It depends on s' and we therefore replace L(8q) by L(s', 8q) in our notation to 
emphasise that it changes as we change the centre s ~ s' of the peaked function f(s, s'). 

If h(s) is sufficiently slowly varying we may therefore approximate the expression for the 
inner product (|3~%j) by 

($|* n > ~ — - - f L L(s',9 )^(s')Ms')ds'. 



Ak 2 cos 60 Jo 

If, however, h(s) is at the same time sufficiently localised around sq, then we can replace this 
in turn by 

4k 2 cos t^o Jo Ak 2 cos & 

This verifies our assertion that (<3?|\l/ n ) and (<p\ip n ) approximately coincide up to a geometrical 
factor 

L(s ,6 l o) 

9 s o,^o) = 77^ j- 

4/c z cos 6^0 

which depends on k and the location of the probe function ip(s) but is independent of if) n (s). 

It should be noted that for scarring around a general periodic orbit we will need to consider 
probe functions which are localised not around a single point (sq, Vq) on the boundary section 
but around a series of points (s t , v t ), t — 0, 1, • • • , T representing the bounces of a periodic orbit 
on the boundary. In this case we consider probe states of the form ip(s) = J2t Ptis) where each 
component is localised in the section around a single point (s t ,Vt). The interior overlap ($|\l/ n ) 
is then a sum of terms of the form given above, each with a geometrical factor g(s t ,v t ) which 
depends on the the angle of incidence of the periodic orbit at the corresponding bounce and the 
length of the orbit segment to the next bounce. We can recover a simple scaling such as given in 
(jH5| if we scale each of the components (pt of ip by the corresponding geometrical factor before 
computing the overlap (or alternatively by redefining the boundary inner product). A detailed 
treatment of this procedure will presumably lead to a formulation in which the boundary 
integral operator is replaced, semiclassically at least, by a unitary evolution (see [21], but do 
not pursue that option here. Instead we note that our explicit calculations have been for a 
two-bounce orbit normally incident on the boundary and in that case the geometrical factor is 
the same for both bounces and no such adjustment is necessary. 
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